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C.  PROBLEM  STUDIED 

Light  scattering  by  spherical  and  nonspherical  particles  and  effect 
of  scattering  on  atnosphet ic  propagation  and  remote  sensing. 


D.  MOST  IMPORTANT  RESULTS 


1.  Perturbation  Approach  to  Light  Scattering  by  Nonspherical  Particles. 

We  have  performed  a  series  of  calculations  applying  the  first 
order  perturbation  theory  to  the  scattering  of  electromagnetic  waves 
by  nonspherical  particles.  The  surface  of  the  nonspherical  particle 
was  parameterized  in  the  form 

r  =  r g  (  1  +  ef  (  6  ,  <J>)  )  , 

where  rg  is  a  radius  of  unperturbed  sphere,  f(0,<t>)  describes  the 
shape  of  nonspherical  particles,  and  the  perturbation  parameters 
e«  1.  We  have  used  several  different  forms  of  the  shape  function 
f(6,4>)»  and  we  have  compared  the  numerical  results  obtained  using  the 
perturbation  approach  with  those  using  the  extended  boundary  condition 
method.  In  this  way  the  region  of  applicability  of  the  first  order 
perturbation  method  was  established. 

2.  Light  Scattering  by  a  Pair  of  Conjugate  Nonspherical  Particles. 

We  define  a  pair  of  conjugate  nonspherical  particles  as  a  pair  of 
particles  that  are  a  mirror  image  of  each  other.  Using  the  first 
order  perturbation  theory  we  have  learned  that  the  sum  of  scattering 
or  extinction  cross  sections  of  these  two  nonspherical  particles 
is  e^ual  to  a  scattering  (or  extinction)  cross  section  of  a  suitably 
defined  sphere. 
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Light  scattering  by  a  pair  of  conjugate  nonspherical  particles 

Petr  Chylek,*  .1.  T.  Kiehl,*  and  A.  Mugnai  * 

Center  for  Earth  and  Planetary  physics.  Harvard  Cnucrsity,  Caml.rulpe.  Massachusetts  0JI3S 
and  National  Center  for  Atmospheric  Research,  boulder,  Colorado  H03U7 
(Received  1  February  1979) 

Applying  the  first -order  perturbation  theory,  we  have  derived  a  theorem  which  states  that,  under 
specific  conditions,  the  sum  of  the  scattering  matrices  of  two  nonspherical  particles  can  be  replaced 
by  the  scattering  matrix  of  a  sphere.  Consequently,  a  polydispersion  of  such  nonspherical  particles 
can  be  replace.!  by  a  polydispersion  of  spheres,  without  changing  scattering  characteristics  of  the 
polydispersion.  This  implies  the  nonuniqueness  of  the  inverse-scattering  problem.  To  verify  the  valid¬ 
ity  of  the  theorem,  the  differential  scattering  cross  sections  of  several  nonspherical  rotationally 
symmetric  particles  have  been  calculated  using  the  extended  boundary  condition  method.  The  re¬ 
sults  show  that  the  theorem  is  satisfied  with  accuracy  expected  from  the  first-order  perturbation 
theory. 


INTRODUCTION 

Light  scattering  by  small  particles  has  found  many  appli¬ 
cations  in  astronomy,  biology,  chemistry,  meteorology, 
physics,  and  other  branches  of  science.  In  some  situations 
the  scattering  particles  are  spherical  (for  example,  fog  and 
cloud  droplets)  and  in  other  situations  they  are  nonspherical 
(dust  particles  in  the  Earth’s  or  the  Marti  ai  atmosphere,  dust 
grains  in  interplanetary  and  intergal.ictic  space,  macromole¬ 
cules,  cells,  etc.). 

Because  the  scattering  problem  ran  be  solved  analytically 
only  for  very  few  highly  symmetric  shapes  and  because  the 
analytical  solution  for  spherical  scatter!  rs  is  much  simpler 
than  for  other  shapes,  in  many  practical  applications  the  Mie 
theory  (scattering  bv  sphercal  particles)  is  used  regardless 
of  the  real  shape  of  the  scat  te.  in, -:  parti  ,es.  Such  a  replace¬ 
ment  of  a  nonspherical  poivdisp.-rsi  :\  of  particles  by  a 
spherical  polydisper.sion  may  he  acceptable  in  some  cases; 
however,  in  other  cases  it  may  lea  I  to  considerable  dis¬ 
crepancies,  depending  on  the  exact  shapes  of  the  particles 
pres  rt.  Unfortunately,  there  are  no  known  criteria  to  decide 
under  which  conditions  and  with  what  aci  uracy  nonspherical 
particles  can  he  replaced  by  .-pit.  rival  particles.  It  is  the 
purpose  of  this  report  to  show  that,  under  well-specified 
conditions,  nonspherical  part  id-  x  can  he  replaced  hv  spherical 
ones  without  significantly  d.  ic  hr.;  any  of  the  scattering 
characteristics  like  scattering,  ah-  i  p;  V.n  or  extinction  cross 
sections,  differential  scattering  .  Motion,  or  the  state  of 
polarization  of  the  scattered  r.idi  uio". 

CONJUGATE  PARTICLES  THEOREM 

Consider  for  simplicity  a  non..;- here  il  particle  that  can  he 
obtained  by  the  deformation  of  i :  phe:  c  t  !■  ig.  1 )  in  such  a  way 
that  the  depth  t  of  the  layer  it.  which  the  deformation  occurs 
is  equal  to 

(  =  '•’/<,  (1) 

where  r  is  the  radius  of  the  o'iyinal  taic’.-forrped  sphere  and 
t  is  a  number  considerably  sm.i  ler  than  one.  We  should  re¬ 
alize  that  even  if 

<  «  1,  (2) 

the  thickness  t  does  not  h  ,vi  to  he  very  small.  For  example, 
fort  =  0.1,  the  thickness  r  is  cqu.i!  to  PO'.'  of  the  sphere’s  ra¬ 


dius.  It  is  convenient  to  parametrize  the  surface  of  such 
nonspherical  particles  by 

r,  =  r[l  +  tf(0,  0)].  (3) 

where  0  and  <f  are  the  usual  spherical  angular  coordinates  and 
the  function  f(0,  if)  specifies  the  deviation  of  the  particles  trom 
a  sphere  in  the  direction  given  by  angles  l)  and  if. 

Let  us  define  what  we  will  call  a  pair  of  nonspherical  par¬ 
ticles  conjugate  with  respect  to  a  sphere  of  radius  r.  Consider 
two  nonspherical  particles:  particle  A  and  particle  B.  If  the 
surface  of  these  two  particles  is  described  by  the  equations 

Ars  =  r[l  +  e/(0,  0)1,  (4) 

ft r,  =  r[l  -  e/(0,  if)},  (5) 

where  Ar,  is  the  magnitude  of  the  radius  vector  of  an  f.rbh  rary 
point  on  the  surface  of  particle  A  and  p,rs  is  the  magnitude  of 
the  radius  vector  of  a  surface  point  on  particle  B,  we  w.,1  say 
that  the  nonspherical  particles  A  and  B  are  conjugate  with 
respect  to  a  sphere  of  radius  r.  We  can  also  say  that  the 
nonspha.ical  particles  A  and  B  are  conjugate  with  respect  to 
a  sphere  of  radius  r  if  particle  B  is  the  mirror  image  of  particle 
A  with  icspect  to  a  spherical  mirror  of  radius  r  and  vice  versa. 
A  few  examples  of  conjugate  nonspherical  particles  are  shown 
in  Fig.  2. 

Tile  electromagnetic  field  E,  scattered  by  a  given  par’icle 
(spheiical  or  nonspherical)  can  he  written  in  the  form 

E,  =  S  E„  ((',) 

where  K,  is  the  electric  vector  of  the  incoming  radiation  and 
S  is  a  2  X  2  scattering  matrix,  the  form  of  which  is  determined 
by  the  properties  of  the  scattering  particle,  such  as  the  par¬ 
ticle’s  shape,  size,  and  material.  To  determine  explicitly  the 
form  of  the  elements  of  the  scattering  matrix  S,  one  has  to 
solve  Maxwell’s  equations  for  the  boundary  conditions  sat¬ 
isfied  on  the  surface  of  the  scattering  particle:  closed  form 
analytical  solutions  have  been  derived  only  lor  ..j.lierical 
particles,1  spheroids,”  and  i. -.finitely  long  cylinders, ■’  raid  it  is 
doubt lul  that  a  closed  loin  analytical  solution  cat,  exi .:  for 
Rn  arbitrarily  shaped  particle  as  considered  in  this  study. 

Since  the  parameter  f  appearing  in  Eqs.  (3)  ~(,M  i>  consid¬ 
erably  smaller  than  one,  however,  first-order  perturbation 
theory  ran  he  applied  to  the  considered  case  of  light  .ering 
by  nonspherical  p.ntii  les.  It  follows1-'’  that  if  light  s. ;  t-ring 
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FIG  t.  Consideration  of  a  nonspherical  particle  whose  deviation  from 
a  sphere  of  radius  r  occurs  only  in  a  spherical  shell  ol  thickness  t  with  inner 
and  outer  radius  equal  to  t\  1  -  t)  and  4  1  +  t).  respectively. 


by  a  spherical  particle  is  described  by  the  scattering  matrix 
S,  then  the  scattering  by  a  nonspherical  particle  .4,  as  specified 
by  Kq.  (4),  is  described  by  the  scattering  matrix 

„S«S  +  <S,,  (7) 

where  the  matrix  S|  is  independent  of  e.  Since  the  difference 
in  the  mathematical  description  of  nonspherical  particles 
conjugate  to  each  other  with  respect  to  a  sphere  of  radius  r  is 
only  the  sign  of  the  parameter  e  [compare  Eqs.  (-1)  and  (5)1,  the 
scattering  matrix  hS  describing  light  scattering  by  a  particle 


B  (which  is  conjugate  to  4  with  respect  to  a  sphere  of  radius 
r)  is  given  by 

„S  =  S  -  cS,.  (8) 

By  adding  Eqs.  (7)  and  (8),  we  obtain  a  basic  theorem — 
which  we  will  call  the  conjugate  particles  theorem — allowing 
us  to  replace  conjugate  nonspherical  particles  by  spherical 
particles  without  changing  the  observable  scattering  charac¬ 
teristics  of  the  system.  We  obtain 

^S  +  ,;S  =  2S,  (9) 

which  says  that  the  sum  of  the  scattering  matrices  of  two 
nonspherical  particles  4  and  ft,  conjugate  with  respect  to  a 
sphere  of  radius  r,  is  equal  to  twice  the  scattering  matrix  of  a 
sphere  of  radius  r. 

Since  c  «  1,  the  same  additive  property  is  valid  for  experi¬ 
mentally  observable  quantities  like  the  differential  scattering 
cross  section  (phase  function),  absorption,  scattering  and 
extinction  cross  sections,  and  polarization  (we  neglect  terms 
proportional  to  f2  and  higher  powers  of  f ). 

Remember  that  the  sum  ,*S  +  «S  in  Eq.  (9)  means  the  sum 
of  two  single-particle  scattering  matrices  (not  a  two-particle 
scattering  matrix).  Also,  we  consider  the  case  of  incoherent 
scattering  where  distances  between  individual  particles  art- 
much  larger  than  the  wavelength  of  the  scattered  radiation. 
Therefore,  the  sum  is  not  a  function  of  the  distance 

between  particles  4  and  Ft. 


FIG.  2.  Several  examples  ol 
conjugate  nonspherical  particles. 
Nonspherical  particles  A  ("Olid 
line)  and  B  (dashed  line)  a-c-  con¬ 
jugate  with  respect  to  a  sr'  ■  vo  ol 
radius  r.  Particles  show:-  in  (  i) 
and  (b)  are  rotationallv  symmetric 
and  resemble  prolate  and  va.'ale 
spheroids  very  closely.  To  a 
given  nonspherical  particle  A.  ti  e 
conjugate  particle  can  be  direr  a 
particle  6  or  a  particle  9*  de¬ 
pending  on  the  radius  ol  the 
sphere  with  respect  to  wtieh  the 
particles  are  conjugate  [(c).  (J)j. 
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FIG.  3.  Equation  (11)  satisfied 
with  an  error  2  [Eq.  (IT)]  for 
considered  four  pairs  ol  conjugate 
nonspherical  particles  in  Eqs 
(lOaMlOd).  For  <  =  0  01  the 
error  in  all  cases  is  smaller  than 
0.1%.  For  f  =  0. 1,  the  error  is 
smaller  than  10%.  In  cases  (a) 
and  (b)  the  error  i  chafes  the 
sign.  In  all  cases  an  absolute 
value  of  6  is  plotted  and  the  cor¬ 
responding  parts  ol  the  cur.es  are 
denoted  by  the  (+)  and  (— )  signs. 
The  solid  curves  show  the  error 
tor  the  case  when  particles  are 
oriented  in  such  a  way  that  their 
axis  of  rotation  is  parallel  to  the 
direction  ot  the  incoming  light,  and 
the  dashed  curves  show  n  for  ihe 
case  of  perpenrfcular  or  env-'in 
(axis  of  rotation  perpend. ouh'ir  to 
the  direction  of  die  incorr  Irg  I  ght 
and  perpendicular  to  the  scatter¬ 
ing  plane  chosen  to  be  the  x-z 
plane).  The  incoming  radiation  is 
propagating  in  the  direction  ot  the 
positive  z  axis  and  is  linearly 
polarized  with  theelectri;  vector 
making  an  angleof  45“  with  the  x 
axis. 


NUMERICAL  RESULTS 


To  verify  the  validity  of  the  derived  theorem,  we  have  cal¬ 
culated  the  differential  scattering  cross  sections  of  several 
nonspherical  axially  symmetric  particles  using  the  extended 
boundary  condition  method.6-7  We  considered  four  sets  of 
nonspherical  conjugate  particles  given  by  the  equations 

r,  =  r[l  ±  (  sin20]. 

(10a) 

r,  =  r[l  ±  «  T2  (cos'))], 

(10b) 

r,  =  r[l  ±  (  T ,  (cos#)]. 

(10c) 

r,  **  r[l  ±  (  .",(cos  0)], 

(lOd) 

where  T„(c os  0)  is  the  n  th  Cbebvshev  polyno 

tnial  of  the  first 

kind,  and  where  Pi(cosfl)  is  the  fourth-order  Legcrdre  poly¬ 
nomial.  While  (10a)  and  ( 10b)  are  convex,  particles  described 
by  F,qs.  (10c)  and  (l(hl)  are  concave.  Cross  sectional  views 
of  the  particles  (10a)  and  (10b'  ;  re  schematically  shown  oil 
Figs.  2(a)  and  2th).  Let  us  con  uh-r  the  relation  between  the 

|5fi2  J.  Opt.  Soc  Am.,  Vol.  fit), 11,  November  197;) 


differential  scattering  cross  sections  of  the  foim 

Ao  +  no  =  2ro  (11) 

following  from  the  theorem  [Kq.  (9)]  for  the  caser  «  1  (in  the 
above  notation  rn  stands  for  the  differential  scattering  cross 
section  c!  a  sphere  of  radius  r).  In  Fig.  3  the  percentage 

error 

,  2,(7  —  (  j  (7  +  no) 

6  =  — - — - —  X  100%  (12) 

A(T  +  „(7. 

is  plotted  a.;  a  function  of  the  scattering  angle  for  two  duY-.-rent 
orientations  of  the  scattering  particles.  In  general,  we  can  sav 
that  for  (  =  0.01  (deviation  from  a  sphere  occurs  in  an  outer 
layer  of  thickness  equal  to  2%  of  the  radius),  the  theor  .m  is 
satisfied  with  an  accuracy  h  tier  than  0.1%.  With  (  =  0,0'.  (t 
equal  to  in%,  of  the  radiiLs),  the  accuracy  is  than  2%;  and  finally 
with  e  ~  J.I0  (t  equal  to  20V  of  the  radius),  the  amir  icy  is 
alwa  s  better  than  10%  .  This  accuracy  is  in  agreement  with 
what  nr. i-  should  expect  from  neglecting  terms  pro)  or'i.mal 
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to  t-  with  respect  to  i.  We  conclude  th  it '  ::e  numerical  results 
confirm  the  validity  of  the  conjugate  particles  tlieorem. 

APPLICATIONS 

The  derived  theorem  has  several  practical  consequences. 
Whenever  we  have  a  pulyd:-[>‘-rsicm  of  nonspherical  particles 
such  that  for  each  nonspherical  particle  of  the  polydispersion 
there  exists  another  nonspherical  p  rticle  in  the  same  poly¬ 
dispersion  that  is  conjugate  l  "  the  pr  vious  one  with  respect 
to  a  sphere  of  radius  r,  then  the  con-  ■  ierc  ri  polydispersion  of 
nonspherical  particles  can  he  repine  ;  hy  a  polydispersion  of 
spheres  (each  pair  of  nonspherical  particles  conjugate  with 
rcspoct  to  a  sphere  of  radius  r  is  replan  d — according  to  the 
theorem — by  two  spheres  of  radius  r)  without  a  change  of 
scattering  characteristics  of  the  polydi-nersion.  It  means  that 
within  the  accuracy  of  order  <  all  scattering  characteristics  of 
both  polydispersions  are  identical  and  there  is  no  experi¬ 
mental  measurement  that  can  distinguish  between  these  two 
polydispersions  on  the  has  -  of  their  scattering  properties. 
This  shows  directly  the  nontiniquet;;ss  of  the  inverse-scat¬ 
tering  problem.  From  the  scattering  characteristics  measured 
with  a  definite  experimental  error  cannot  deduce  unam¬ 
biguously  the  size  distribution  of  the  .-mattering  particles  (even 
if  we  assume  that  the  index  of  rt  frac  tion  is  known).  Even 
more,  we  are  not  able  to  determine  whether  particles  are 
spherical  or  nonspherical.  Of  course,  with  increasing  exper¬ 
imental  accuracy  the  range  of  ambiguity  decreases. 

If  a  nonspherical  particle  is  self-conjugate  (conjugate  to 
itself  in  a  different  orient etionl  with  respect  to  a  sphere  of 
radius  r,  then  the  scattering  mat.-,  .:  of  this  particle  averaged 
over  all  orientations  is  c.nnl  o  the  scattering  matrix  of  a 
sphere  of  radius  r. 

If  wc  write  the  conjugate  'tick,  theorem  of  Eq.  (9)  in  the 

form 

/iS  =  2  S  -  ,,S,  (13) 

we  can  see  that  the  scattering  itr\  of  the  nonspherical 
particle  H  can  be  calculated  ••me  v.->  know  the  scattering 
matrix  of  the  nonspheF-.tl  par:  ..  ie  .-1  conjugate  to  B  with  re¬ 
spect  to  a  sphere  of  radius/-.  let  u-as-ume  that  the  scattering 
matrix  of  the  nonspherical  p.ir’a  1  •  -.•>•>,■  n.  We  can  make, 

for  example,  a  model  of  the  ..  .inn  ie  .1  and  perform  micro- 
wave  scattering  experiment  '•  cu  t’  mono  its  scattering 


matrix.  Then  we  can  choose  a  set  of  spheres  with  radii  rh 
r  , ,  r.|,  .  .  .  ,  r,  and  determine  the  shapes  of  the  nonspherical 
particles  /tj.  /E, conjugate  to  the  nonspherical  particle 
A  with  respect  to  the  spheres  of  radii  r  t,  /••.,  r ... Tiien 
the  scattering  in.ilrjcta  for  the  nonspherical  particles  lly, 

H  i,  /?-, . H,  can  be  calculated  from  Eq.  (13).  Thus,  one 

experimental  measurement  can  provide  scattering  matrices 
for  several  differently  sh.aped  nonspherical  particles. 

Finally,  one  can  similarly  define  nonspherical  particles 
conjugate  with  respect  to  a  surface  other  than  a  sphere  'fur 
example,  with  respect  to  an  oblate  or  a  prolate  spheroid  i  and 
then  replace  a  polydispersion  of  irregular  conjugate  particle- 
by  a  polydispersion  of  spheroidal  particles  the  scattering 
property  of  w  hich  can  be  calculated  analytically. 
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ABSTRACT 

A  linear  relationship,  in  v -enJent  of  the  form  of  the  size-distribution,  between  extinction  at  wave¬ 
lengths  around  A  -  lljtr.i.  •  'priori  around  A  =  3.8  and  9.5  Mm.  and  liquid  water  content  of  atmospheric 
fogs  has  been  verified  u-rr-c  341  droplet  size  distribution  measurements  made  under  a  variety  of 
meteorological  conditions.  I  he  results  suggest  that  integrated  liquid  water  content  along  a  path  in  fog 
can  be  determined  from  p-vsuremont  of  CO.  laser  (A  =  10.6  ^m)  transmission  along  the  path,  and  that 
liquid  water  content  at  a  t-  o-  cular  point  in  fog  can  be  inferred  from  in  situ  measurement  of  fog-droplet 
absorption  with  ?  deute:  am  fluoride  la«cr  (A  =  3.8  gm)  or  a  suitably  tuned  CO.  laser  (A  =»  9.5  pm) 
spectrophone. 


1.  Introduction 


It  has  been  recently  shown  (Chylck,  1978)  that  a 
linear  relationship,  independent  of  the  form  of  the 
size  distribution,  should  exist  between  the  infrared 
extinction  around  A  =  11  »-.m  and  the  liquid  water 
content  of  toys.  The  relation  can  be  written  in  the 
form 


I7?-  H\ 

2  f>K 


(1) 


where  cr,.  is  the  volume  extinction  coefficient  meas¬ 
ured  at  the  wavelength  \,  IV  th..  water  content. 
(>  the  density  of  water,  and  the  coefficient  c  is  equal 
to  the  slope  of  a  straight  line  that  approximates 
the  extinction  efficiency  curve  (J,  m  ,\)  by 

QA X.k)-  c(\h,  (2) 


where  the  size  parameter. v  is  defined  by  the  ratio  of  the 
particle  circumference  to  the  v.  .-.sc length.  An  ap¬ 
proximate  value  of  the  coclnciciit  <  (A)  a!  a  -  II  n  m 
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isc  0.31 .  T.te  conditions  under  which  the  approxi¬ 
mation  (2)  arc  valid  and  the  derivation  of  the  rela¬ 
tion  (I)  have  been  discussed  elsewhere  (.Chylek, 
1978). 

In  this  paper  we  verify  the  validity  of  relation  (I) 
by  calculating  the  volume  extinction  coefficient  cr,. 
and  the  liquid  water  content  W  for  341  different  fog 
droplet  size  distributions  (Garland,  I97T,  Kumai, 
1973;  Garland  ot  ul..  1973;  Kunkcl,  1971:  Roach 
et  ul.,  1976;  Pinnick  it  ul.,  1978)  measured  under 
various  meteorological  situations. 

We  also  show  that  a  linear  relationship,  similar 
to  (1),  Csisis  between  the  infrared  absorption 
coefficient  in  the  spectral  teginie  \  =  3. 5-5. 3  jcm. 
8-10  /am.  and  lire  liquid  water  content  of  fogs. 
Thus  for  example,  the  absorption  coefficient  of  fogs 
at  X  =  3.8  jun  is  uniquely  related  to  their  extinction 
and  absorption  at  X  -  10  /am 

2.  Selected  f(  g  size  distributions 

Relatively  few  reliable  measurements  of  fog 
droplet  size  d-'stiibutions  have  been  made,  particu- 
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larly  for  which  numerical  data  are  available.  The 
fog  measurements  used  here  we  judge  to  be  reliable 
and  were  chosen  to  represent  a  wide  range  of  fog 
conditions  ranging  from  maritime  and  continental 
advection  fogs  [Kumai,  1973;  km  ..cl,  197 1;  and  part 
of  Garland’s  (1971)  work]  to  ininml  radiation  fogs 
(Garland.  1971;  Garland  ct  at..  1975;  Roach  el  til., 
1976;  Pinnick  et  al.,  1978).  The  early  work  on 
evolving  fogs  near  the  Atlantic  Ocean  in  France 
and  stable  inland  fog  and  haze  near  Paris  by  Arnulf 
et  al.  (1957)  was  not  used  because  it  was  not 
possible  to  obtain  true  droplet  distributions  from 
their  figures.  Arnulf  et  a!.  (1957)  captured  droplets 
on  spider  threads  (for  which  the  capture  coefficient 
depends  on  droplet  size)  but  gives  only  the  uncor¬ 
rected  droplet  distribution  data.  Results  of  the 
pioneering  work  of  May  (1961)  using  a  specially 
designed  two-stage  impactor  were  not  used  since 
the  numerical  data  are  no  longer  available  (May, 
private  communication,  1978).  Measurements  of 
valley  fog  drop  sizes  obtained  by  exposing  gelatin- 
coated  slides  to  a  stream  of  foggy  a.ir  by  Pilie  et  al. 
(1975)  were  deemed  not  eic-ible  sit.ee  the  distribu¬ 
tions  were  normalized  to  simultaneous  measure¬ 
ments  of  extinction  coefficient  derived  by  a  trans- 
missometer.  Fog  drop  measurements  made  with  a 
light-scattering  counter  by  Fldrid  /-•  (1961)  were  not 
utilized  since  we  suspect  errors  in  his  measure¬ 
ments  due  to  non-isokinetie  sampling.  The  inlet  of 
his  counter  was  only  I  cm  in  diameter  and  became 
wet  during  the  sampling  ;  races-.  In  addition,  some 
unexplained  differences  '.ere  cait-ed  by  a  dilution 
apparatus  Mldridge  used  !  U  high  droplet  concentra¬ 
tion  conditions. 

Three  different  sampling  techniques  were  em¬ 
ployed  to  obtain  the  fog  size  distributions  utilized 
in  this  study;  impaction,  holographic  and  light  scat¬ 


tering. 

Garland  (1971),  Garlan  I  et  id.  (1973)  and  Roach 
ct  al.  (1976)  used  a  nv  Jified  two-stage  Caselia 
impactor  designed  by  May  (1961)  mounted  hori¬ 
zontally  in  a  wind  tunnei  to  provide  isokinetic 
sampling.  Corrections  to  the  collection  efficiency 
based  on  the  penetration  cm  so-  of  the  impactor 
(May,  1945)  were  applied  n>  ' 
vice  is  sensitive  to  dropk.s  o 

Kumai  (1973)  also  used  : 
to  measure  advection  fog  dro 
Arctic  Ocean  at  Point  !>..  ■. 
with  a  gelatin-coated  plu 
whose  collection  efficiency  v. 
lion  of  wind  velocity.  Also, 
tion  factors  were  applied  ;c 
combined  methods  yielded 
from  2.2-64  /mi  radius. 

In  general,  the  primary  d 
tion  technique,  besides  i-qt. 
duction,  is  the  uncert  i  uy 
tion  determination  foi  ueu 


:.iw  data.  This  de- 
■  -72  /un  radius, 
o -stage  impactor 
!  :  nied  over  the 
Alaska,  together 
collection  plate 
a  liable  as  a  func- 
droplet  correc- 
n ue  sizes.  T  he 
concentration 

of  the  impac- 
irious  data  re- 
er  concent  ra¬ 
cier  ami  sub- 


micrometer  size  droplets,  in  addition,  the  size  limit 
of  detectability  is  about  0.3  pm  radius  using  con¬ 
ventional  microscopy  techniques. 

A  laser  hologram  technique  was  employed  by 
Kunkel  (1971)  to  measure  droplets  in  advection  fog 
propagating  inland  during  nighttime  at  Otis  /'.ir 
Force  Base,  Massachusetts.  The  hologram  camera 
was  capable  of  sampling  volumes  of  4.5  cm'  at  a 
rate  of  five  samples  per  minute,  in  a  near-isokinetic 
fashion,  and  with  minima!  disturbance  to  the  drop¬ 
lets.  Droplets  with  radius  2-40  pm  were  detected. 
Because  of  the  small  number  of  droplets  (normally 
<50)  in  the  distributions  reported  by  Kunkel,  we 
have  averaged  all  17  of  the  reported  droplet 
distributions  together,. which  in  any  case  repa  vents 
only  a  3-min  interval,  to  obtain  a  single  distribution. 

Finally,  we  have  used  measurements  by  Pinnick 
et  al.  (197S)  of  radiation  fog  and  haze  made  during 
wintertime  in  West  Germany  w'ith  a  commercially 
available  light-scattering  counter  (the  “Knollcn- 
berg”  Class'cal  Scattering  Aerosol  Spectrometer 
manufactured  by  Particle  Measurement  Systems, 
Boulder,  Colorado).  This  device  works  on  the 
principle  that  as  aerosol  flows  through  an  illuminated 
volume,  light  scattered  by  single  droplets  into  a 
particular  solid  angle  is  measured  and  used  to 
determine  particle  size  by  making  pulse-height 
analytes  of  the  lesponse  pulses.  Determination  of 
droplet  size  from  the  response  is  indirect  because  of 
the  dependence  of  the  response  on  factors  other 
than  dropiet  size,  viz.,  droplet  refractive  index  and 
the  lens  geometry  of  the  optical  system.  Particular 
attention  was  given  to  the  calibration  of  this  insO  i  ■ 
ment  using  monodispersc  particles  of  different  s.r  ; 
and  refractive  index.  The  manufacturer's  adv.  - 
tised  calibration  was  not  used.  Rather,  droplet  sme 
distributions  were  determined  by  redefining  the  size 
pulse-height  channels  as  described  in  detail  by  Fin- 
nick  et  al.  (1978).  This  counter  is  sensitive  to  water 
droplets  with  radii  0.23-  16  gun. 

Altogether.  34!  different  size  distributions  were 
used  to  check  the  validity  of  F.q.  (1)  in  the  atmos¬ 
pheric  window  around  A  -  11  /am;  25  fog  distribu¬ 
tions  were  taken  from  Garland  (1971)  and  Roach 
ct  ( d .  (19/9),  6  from  Garland  ct  al.  (1973),  2f'  Horn 
Kumai  (1973),  I  from  Kunkel  (1971),  and  289  tog 
and  haze  distributions  from  Pinnick  ct  al.  (1978). 


3.  Numerical  results  for  log  extinction 

Using  a  Mie  scattering  program  and  index  of  re¬ 
fraction  o'  water  ns  given  by  Hale  and  Querry  ( 1973). 
we  have  calculated  the  volume  extinction  cocflieu  ut 
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with  r  >  14  /rni;  however,  contribution 

generally  did  not  dominate  eiti-  ■:  ilie  extinction 
rr,  or  the  liquid  water  content  IV.  Consequently, 
we  expect  that  the  linear  relationship  (1)  between 
the  volume  extinction  coefficient  rr,  and  the  liquid 
water  content  U'  will  be  reasonably  well  satisfied 
at  A  ---  11  /on.  Results  of  nu  nerical  calculations 
confirming  the  validity  of  Fq.  (1)  at  X  —  11 
arc  shown  in  Fig.  Id. 

On  the  other  hand  at  X  =  0.55  pm  the  approxima¬ 
tion  Q,.  ■---  r.r  is  valid  only  for  water  droplets  with 
r  0.5  pm  (Chyleh.  19/5  .Since  most  haze  and 
all  fog  droplet  size  distributions  me  dominated  by 
droplets  with  r  >  0.5  /.  in,  the  Q.  -  t.v  approxima¬ 
tion  is  not  valid  in  this  case  and,  consequently, 
no  size-distribution-inde:  endent  (elation  between 
the  extinction  and  liquid  w.tt-e:  content  should  exist 
at  X  =  0.55  pm.  The  numerical  results  based  on  the 
measured  distributions  commo  .-.1  to  the  F.q.  (1) 
approximation  in  Fig.  la  c  inis  conclusion. 

Further,  the  approximation  grossly  overestimates 
extinction  for  most  of  the  distributions. 

As  we  change  the  wavelength  bom  X  -  0.55  /am 
to  longer  wavcL np  '  ti  e  O,  ••••  .  npioximation  is 
satisfied  for  larger  droplet  ra>  .  Consequently, 
with  increasing  wavelength  the  i.  don  given  by  Fq. 
(1)  is  becoming  a  more  realistic  approximation  for 
hazes  and  for  fogs.  This  trend  c:u> !  .•  s  .  n  from  Figs, 
lb-lc,  showing  the  numerical  •.••!;t;,.ttious  and  the 
Eq  (1)  approximation,  a.  X  -  and  -I  /an.  Wc 
note  that  with  mere:; s lug  X  the  ’elation  (I)  better 
approximates  the  e.xac'  :>e metical  results  and  finally 
at  X  =  11  pm  (Fig.  la!  the  ••  nlume  extinction 
coefficient  becomes  i  *•.. . ;  endent the  si/e  distribu¬ 
tion  /if/-),  and  the  relation  ( ! )  agteo  s  within  a  factor  2 
with  the  numerical  results  for  ..Is  distributions. 

Closer  inspection  of  the  test:!:  in  Fig.  Id  shows 
noticeably  better  agi  content  :  etween  the  Mie 
numerical  results  and  the  ivh.tioi:  (I)  for  the  size 
distributions  of  Finniek  rr  til.  (lf|7M.  The  reason 
very  likely  has  to  do  with  t!  "act  that  the  size 
distribution  measurements  v ■■  i  ;  ink  rt  al.  are  only 
for  droplets  with  t  u:  c;>  to  r  -  19  pr.i.  whereas 

the  Ciarhmd  i  197 1 ).  Gai la  id  it and  Roach 
cl  til.  (19761  measurements  ate  to;  droplets  with 
radii  up  to  r  ■  72  /nr.:  'lie  liu:  tai  .I0-.')  measure¬ 
ments  are  for  droplets  with  n  to  /  64  /tin; 

and  the  Kttnkel  (197li  tm-nsur-  a  :  is  for  droplets 
with  i.  dii  up  to  r  40  mu.  I  ,  the  maximum 
radius  condition  for  the  »J,  i  appioximation 
leading  to  relation  (iiis  tor  r.  ’  pm  at  X  -■  11 
/tin,  none  of  the  1’itmi.  k  it  ■  uihmions  can 
strongly  violate  this  condition  panicles  with 

/  7*  16  /tin  were  metsuivd.  ’  the  better 
agreement  of  the  numeii:a!  ies;  r  t.he  I’minek 
rt  nl.  distiibution  .  wi'h  .he  :c!  ’ll  at  X  II 

/tm  [and  also  at  X  -  4  /.m  tsc-c  I  >|  may  be  in 


part  a  consequence  of  their  inability  to  measure 
droplets  with  r  >  16  // m. 

Another  qualification  concerning  the  results  in 
Fig.  1  bears  on  our  assumption  that  all  fog  and  haze 
particles  consist  of  homogeneous  water  droplets 
and  have  complex  refractive  indexes  of  water. 
Haze  particles  in  particular  may  contain  a  signifi¬ 
cant  volume  fraction  of  contaminants  such  as  sea 
salt  or  ammonium  sulfate.  The  crucial  question 
here  is  to  what  degree  is  the  particle  refractive 
index  affected  by  such  contaminants?  We  know 
in  the  case  of  fog  that  its  formation  requires 
atmospheric  relative  humidity  RH  to  be  near  10057, 
We  also  know  that  for  the  haze  data  appearing  in 
Fig.  1  the  relative  humidity  w'as  close  to  100(7 
(Pinnick  rt  uL,  1978).  Hand  (1976)  and  Hand 
and  Bullrich  (1978)  have  studied  the  effect  of  rela¬ 
tive  humidity  variations  on  mean  complex  refrac¬ 
tive  indexes  of  maritime  and  urban  aerosols.  Hand 
(1976)  found  that  for  RH  5  95%  the  real  and 
imaginary  parts  of  the  complex  refractive  index 
and  Him  at  X  -  0.55  /im  have  values  1.33  'r 
;;  1.36,  0  •-  /i1;„  v  0.006.  Fxamination  of  Mie  elti- 
cicncy  factois  (_>,(» i.x)  for  refractive  indexes  in  this 
range  shows  our  assumption  that  m  —  1.33-0/ 
in  Mie  calculations  of  extinction  according  to  (3) 
and  the  0,-  -  r.\  approximation  (1)  is  a  good  one.  At 
wavelengths  >X  -  0.55  /on  we  find  from  Hand 
and  Bullrich’.-.  formulas  for  values  of  Rll  >9557 
that  both  maritime  and  urban  aerosol  mean  refrac¬ 
tive  indexes  are  again  not  markedly  different  from 
those  of  water.  For  example,  at  \  --  1 1  /tin  we  pic- 
dicl  fiont  Hand  and  Bullrich’s  formulas  that  the 
real  and  imaginary  parts  of  the  complex  index  a  t 
1.153  v  -s  1.207,  0.0968  =7  /ii3.  0.109,  com¬ 

pared  to  i.-i  ■  1 . 153 -0.095Xf  for  pure  water.  The 
effect  of  '.hose  refractive  index  variations  in  Mie 
calculations  of  extinction  coefficient  according  to  (3) 
and  (I)  a ie  estimated  to  be  not  more  than  10'7. 
For  fog,  an  additional  argument  can  be  made  to 
support  our  assumption  that  the  particle  refractive 
indexes  can  be  approximated  by  those  of  ptue 
water.  1  he  argument  is  that  the  liquid  mass  con¬ 
tent  of  legs  is  on  the  ol  der  of  0.005  g  in  :  orgteatei . 
and  thus  the  volume  Faction  of  any  contaminant 
in  fog  dioplcts  must  necessarily  be  small  so  that 
the  refractive  indexes  must  be  close  to  those  of  watei 

In  older  to  examine  more  closely  the  fog  lesulv 
in  Fig.  1  in  towns  of  fog  tv  pc  we  have  chosen  to 
restrict  our  attention  to  the  data  of  Garland  (I97h. 
Garland  rt  <il.  (1973)  and  Roach  rt  ill.  (1976).  I  lie 
icasons  are  twofold:  fust,  these  measiuemenls  wete 
made  during  ‘mgs  oceanic-,:  under  distinctly  ddlei 
ent  meteoi ologicol  conditions  Altogether,  37  dii br¬ 
ent  logs  v.eic  measined  dmmg  ;i  live-) eat  |vi;o.’ 
uiulei  the  gamut  of  meteorological  conditions  fotmo 
in  England.  Second,  mcasi.i emcn'.s  wete  made  loi  a 
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sufficiently  broad  range  of  particle  sizes  (0.3  pm 
<  r  <  72  pm)  that  error*  in  extinction  and  liquid 
water  content  due  to  the  presence  of  larger  and 
smaller  droplets  are  estimated  to  be  small.  Of  the 
37  measured  distributions,  three  were  not  used  due 
to  nonavailability  of  the  raw  data,  three  because  of 
the  presence  of  ice  crystals  in  the  samples,  and 
five  were  not  used  because  fog  type  was  not 
specified. 

We  have  divided  the  Garland  and  Roach  el  al. 
fog  data,  which  already  appear  in  Fig.  1,  into  two 
classes:  radiation  fog  and  advection  fog.  We  have 
been  cautioned  (Garland,  private  communication, 
1978)  that  although  the  radiation  fogs  clearly 
formed  in  situ  by  radiation  cooling,  some  fogs 
classified  as  advection  type  may  have  been  mature 
radiation  fogs  transported  by  the  wind  from  a  distant 
area  of  formation.  In  any  case,  the  data  are  re¬ 
plotted  according  to  this  classification  in  Fig.  2 
(X  =  0.55  pm)  and  Fig.  3  iX  =  1 1  /am).  At  X  =  0.55 
/am  it  is  evident  radiation  fogs  are  generally  more 
effective  scatterers,  and  hence  more  effective  in 
reducing  visibility,  than  advection  fogs  with  the 
same  liquid  water  content.  To  understand  the  rea¬ 
son  for  this  result,  we  have  picked  a  radiation  fog 
(Fig.  2,  solid  circle)  and  an  advection  fog  (Fig.  2, 
solid  square)  measurement  with  about  the  same 
liquid  water  content,  and  have  plotted  their  differ- 
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Fig.  2.  As  in  Fig.  la  exccp!  *ftat  only  the  fog  of  Ci.tr  I  ami 

( 1971).  Garlnmf  ct  al.  ( 1973)  atA)  Roach  t  t  <•/.  i  !9>*i  arc  shown. 
The  extinction  ami  liquid  w:mr  e  x  tents  cjfci.'latevJ  from  the  26 
measured  size  distributions  in:  divided  according  to  radiation 
(circles)  or  advection  ( square  \  g  We  see  from  the  figure  that 
ru'iation  fogs  are  generally  ..ore  effective  sc.vtcrcrs  thin 
advection  fogs  with  the  same  liquid  water  content  because  they 
contain  more  dioplets  in  the  Vie  resonance  region  that  con¬ 
tribute  a  significant  part  id  ;he  extinct. on  hut  contribute 
only  a  marginal  amount  to  the  iqmd  water  content 


Fig.  3.  As  in  Fig.  2  except  for  K  -  II  ,um.  The  predicted 
relation  between  extinction  and  liquid  water  content  asen  by 
Eq.  (1)  is  shown  by  the  straight  )ine.  Because  the  (7,  "  e\  ap¬ 
proximation  is  generally  better  satisfied  for  radiation  fogs, 
those  points  fall  closer  to  the  straight  line  prediction. 

ential  extinction  coefficients  at  X  =  0.55  pm  vs 
particle  radius  in  Fig.  4.  Because  the  plots  are  m;'  .le 
on  a  linear  scales,  the  areas  under  the  curves  are  a 
measure  of  the  corresponding  extinction  coefficients. 
Thus  the  radiation  fog  extinction  coefficient  (.’5 
km'1)  is  more  than  twice  the  advection  fog  extinc¬ 
tion  coefficient  (15.8  km  ').  For  the  radiation  fog. 
small  droplets,  say.  r  3  gm,  arc  very  numerous  and 
contribute  60 Vr  of  the  extinction  without  making  a 
significant  contribution  to  liquid  water  content.  On 
the  other  hand,  for  advection  fog.  a  broader  size* 
distribution  is  found  and  these  smaller  particles 
contribute  only  6 Cf  of  the  extinction  at  X  =-  0.5f 
pm.  Both  tin*  radiation  and  the  advection  fog  si.’c- 
distributions  strongly  violate  the  maximum  radius 
condition  allowed  in  the  Q,  -  c.v  approximation 
( rm  ■---  0.5  /cm  at  X  -  0.55  pm),  so  there  is  no  reason 
to  expect  a  unique  relation  between  extinction  at 
X  -  0.55  pm  and  liquid  water  content. 

At  X-  II  pm  (Fig.  3),  although  the  Qr~cx 
approximation  is  within  a  factor  of  about  2  for  both 
radiation  and  advection  fog  results,  the  approxima¬ 
tion  is  generally  in  better  agreement  with  the  radia¬ 
tion  fog  results.  The  explanation  is  that  the  (J, 
approximation  is  only  strictly  valid  providing  fog 
droplets  h we r  «-  l4pm(Chylok,  1978)  and  radiation 
fogs  better  satisfy  this  condition  than  do  advection 
fogs.  The  oegree  to  which  this  maximum  '.ulius 
condition  is  violated  can  he  aeterniined  for  two  fog 
examples  from  Fig.  5.  Shown  is  the  differential 
extinction  coefficient  at  X  -  11  pm  for  a  ladietin-. 
fog  (Fig.  3,  ...lid  circle)  and  an  advection  fog  (lag.  *. 
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Fig.  4.  Differential  extinction  coefficients  at  X  =  0.55  /xm  versus  droplet 
radius  for  a  radiation  fog  (solid  circle  in  Fig.  3)  and  an  advection  fog  (solid 
square  in  Fig.  3)  with  about  the  same  liquid  water  content.  The  areas  under  the 
curve*  are  a  i  .easure  of  the  total  extinction  coefficients  (35  km'1  for  the 
irJiat  tm  f.-.g  versus  15.8  km'1  for  the  advection  fog).  For  the  radiation  fog. 
('! of  tiit-  extinction  is  contributed  by  droplets  having  r  3  gm,  compared 
to  only  (■(.  t  r  the  advection  fog. 


solid  square)  again  1  .vita:  about  the  same  liquid 
water  content.  For  the  t;  !i;  lion  fog  69%  of  the 
extinction  arises  from  divaf-is  with  radii  less  than 
the  maximum  value  r,„  14  p in,  while  for  the 

advection  fog  this  value  drops  to  33%.  A  survey  of 
all  the  differential  cxtiiic-  on  coefficient  versus 
particle  radius  graphs  sir.41  .•  those  show  n  in  Fig.  5 
for  the  radiation  fogs  c.  fiat!  u!  (1971).  Garland 
el  dl.  (1973)  and  Roach  .  t  ( 19/6 I  in  Fig.  3  show 
that  in  all  cases  extinction  .  •  A  •  1  inn  is  dominated 
by  droplets  with  r  -  '  14  /-:.:.  <");.  the  other  hand,  a 
survey  of  the  differential  extinction  coefficient 
graphs  for  the  advection  i-  ;  ; csenied  in  big.  3, 
and  also  the  advection  *'  usnai  (1973),  show 

that  for  ttboul  one-h.di  !i  ibutions.  extinc¬ 
tion  at  A  --  11  ytm  is  ■  ’  by  droplets  with 

r  >  14  /on.  Since  •!..  .  p  efficiency  factor 

for  these  larger  path  restimated  by  the 

Q,  c.i  iipproximatiiu.  .  :  MuCiiylek,  1978). 

the  numerical  calculation  of  .  .tsnetion  lot  these 
advection  fogs  I, ill  l>c!o\  !:  *  '■  prediction  in 

Fig..  3.  l  or  the  remaining  li  .ivcction  fogs, 

droplets  wit  hr  •  14 /v  ■  den  i  action  and  the 

points  fall  within  20%  if  i.  I)  piediction. 


Thus,  while  droplets  with  r  =s  14  gm  dominate 
extinction  at  X  —  11  /im  for  radiation  fog,  this  is 
not  always  the  cace  for  advection  fog,  where  the 
presence  of  larger  droplets  partially  destroys  the 
size  distribution  independent  linear  relation  (1). 

4.  Q„  f'x  approximation  for  absorption 

Realizing  that  the  (J,  ■-  c.v  approximtition  wotkv 
reasonably  well  for  fog  at  A  11  /am,  and  for  haze 
at  shorter  wavelengths,  we  checked  to  see  if  a  simi¬ 
lar  approximation  for  fog  droplet  absorption  midst 
hold  in  the  A  -  3-5  /an  and  A  =  8-12  /am  atmos¬ 
pheric  window  spectral  regions. 

The  absorption  coefficient  «>■„  for  ;i  polydispersiou 
of  droplets  described  by  the  size  distribution  n(/  ) 
is  given  by 

a, ,  =  r,r(J„(  in ,  v)n(r)<lr,  (3) 

where  (J„i*w.d  is  the  Mie  efficiency  factor  Ini 
absoiption  im  a  water  dioplet  with  refractive  indc' 
»i(A)  and  size  parameter  a  2~r!\.  Plots  of  th  - 
efficiency  facioi  for  absorption  Q„  vs.t  for  A  3.8, 
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Fig.  5.  Differential  extinction  coefficients  at  X  =  II  fim  versus  droplet  radius  for  a 
radiation  fog  (solid  circle  in  Fig.  3)  and  an  advection  fog  (solid  square  in  Fig.  3) 
with  about  the  same  liquid  water  content.  The  total  extinction  for  the  radiation  fog  is 
22.6  km-1  compared  to  15.0  km"'  for  the  advection  fog-  The  fraction  of  extinction 
contribut'd  1  y  droplets  with  r  s  14  #rm  (the  maximum  value  allowed  in  the  Q,  -  cx 
approximation)  is  69Cc  for  the  taxation  fog,  decreasing  to  33 Cl-  for  the  advection  fog. 
Thus  the  prediction  between  extinction  (at  X  -  II  pirn)  and  liquid  water  content 
(given  by  tq.  (1)  and  shown  in  Fig.  3)  is  a  better  approximation  for  radiation 
fogs  than  for  advection  fogs. 


9.5  pm  are  shown  in  Figs.  6  and  7.  Again  we  find, 
as  Chylck  ((978)  found  for  extinction,  that  Qa  can 
be  well  approximated  by  (70(.t,A)  =  provid¬ 

ing  a:  xm. 

Using  this  linear  approximation  for(7„  in  (5)  gives 


f  “irr1 

J 


O 

oc.  0.5 
o 
»- 
o 

<  0.4 


Thus,  explicit  dependence  on  the  size  distribution 
disappears  and  leads  to  the  absorption  coefficient 
being  linearly  related  to  liquid  water  content  IF 
according  to 
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As  in  the  case  of  cxtinc'ion  the  restriction  that 
x  «  x m  need  not  be  strictly  satisfied,  but  water 
droplets  with  radii  greater  than  the  value  r,„ 
-  Axm/2rr  must  not  contribute  excessively  to  either 
absorption  or  liqui-l  a.uu  content.  Numerical 
values  of  the  maximum  ra  lii  r,„  at  various  wave¬ 
lengths  \  as  well  as  the  s'ope  < '  of  a  sti  night  line 
approximating  Q„  for  \  a,,,,  are  given  in  lablc  I. 


Fic. .  6  The  cipciency  factor  for  absorption  (J,,  for  water  versus 
droplet  si/e  parameter  i  at  a  waw length  A  -  _V8  /mi  (index 
of  ictr.iction  MW-II.IHI.U'i  The  efficiency  factor  c.m  b 
approximated  '>  a  sti.ught  line  <’i  providing  »  ?  r.s. 
lire  approxim.itior  overestimates  the  exact  value  of  (>„  I  r 
some  si/e  parameters,  hut  underestimates  it  for  other  v 
These  two  firms  tend  to  cancel  leading  to  the  absorption 
coefficient  hemp  linearly  related  to  liquid  water  content  ac 
cording  to  eqiia  t -in  (7). 
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Fio.  7.  As  in  Fig.  6  except  for  X  -  9.5  pm  (index  of  refraction 
of  water  m  =  1.243-0.0443/).  The  efficiency  factor  for  absorp¬ 
tion  can  again  be  approximated  by  a  straight  line  Q„  =  c'x 
providing  x  s  i., 

Also  given  are  values  of  the  quantity  3-r '/2Ap  which 
if  multiplied  by  the  liquid  water  content  give  the 
absorption  coefficient.  In  order  that  these  predicted 
relationships  (and  maximum  radius  conditions)  be¬ 
tween  droplet  absorption  am!  liquid  water  content 
may  be  compared  to  the  corresponding  results 
for  extinction,  the  values  of  t„,  r,„ .  etc.,  for  the 
Qr  =  c.r  approximation  (  ■  > y ! s' k .  ll)78)  for  extinc¬ 
tion  also  appear  ir  Table  I.  The  values  of.v,„,  r,„ 
and  c  for  extinction  ue  sometimes  slightly  differ¬ 
ent  from  those  of  Clef  Ink  !  Iv  ;•)  since  there  is  some 
leeway  in  the  subject  >\ ■:  prve.lurc  fo,  approximat¬ 
ing  the  efficiency  factor  tj.t.v)  by  a  linear  function 
of  size  parameter  v.  We  c.  from  the  table  that  the 
limiting  radius  depend,  strongly  on  wavelength 
and  is  in  general  marke  ils  ciitTcient  for  absorption 
and  extinction  at  a  particular  wavelength.  An  ex¬ 
ception  is  at  A  -  9.5  //in,  where  the  limiting  radii 
are  rm  -  13  pm  for  absorption  as  compared  to  r,„ 
-  12.5  pm  for  cxtinctimi.  Tims  for  fog  droplet 
distributions  that  have  radii  r  <  12.5  gin  (i.e.,  most 
radiation  fogs)  we  can  predict  from  Table  1  that  ab¬ 
sorption  contributes  2lKc  of  thee-  (inctionat  A  =  9.5 
/tin,  independent  of  the  foim  of  the  droplet  size 
distribution.  However,  there  is  obviously  no  unique 
relation  between  log  uosoipiion  and  extinction  for 
all  wavelengths. 

Since  we  have  verified  the  (_\  --  r.v  approximation 
for  extinction  at  X  =•  !i  pm.  which  icquires  that 
droplets  have  maximum  i.aii  is  -  14  ,im,  is  ade¬ 
quate  for  atmospl.erir  log.  .. c  suspect  the  <J„ 
--  c'x  approximation  fo*-  a'  -  ;  . -non  is  adequate  for 

fog  at  selected  wavelength' a'  . >  havin':  /-,,  -  14  pro. 
From  Table  I  we  see  t!  a:  v. lengths  A  -  3.8,  4. 
5.3  and  9.5  pm  either  sali-f.  i-r  nearlv  satisfy  this 
criterion.  Therefore.  ".c  mi  ■  expect  a  linear  icl.i- 
tion  between  fog  absoioii*-;  -r,d  fog  liquid  water 
content  aeeoiding  to  (7...  m  ;  nlenl  of  the  dioplet 
size  disti ibution,  for  thesi;  iilar  wavelengths. 


To  check  this  contention  ue  have  calculated  the 
volume  absorption  coefficient  </„(A)  using  a  Mie 
scattering  program  according  to  Eq.  (5)  for  the 
previously  mentioned  341  fog  and  haze  size  distribu¬ 
tions  at  several  different  wavelengths.  As  in  the  ex¬ 
tinction  calculations  we  have  assumed  particle  re¬ 
fractive  indexes  of  water,  and  have  thus  ncplesled 
refractive  index  differences  that  might  be  caused 
by  the  presence  of  contaminants  such  as  sea  salt 
and  ammonium  sulfate.  The  results  at  A  --  3.8, 
9.5  pm  plotted  as  a  function  of  fog  liquid  water  con¬ 
tent  together  with  the  Q„  -  c'x  approximation  (7) 
arc  shown  in  Figs.  8  and  9.  We  see  that  although 
the  numerical  results  are  slightly  better  approxi¬ 
mated  at  A  =  3.8  pm  as  compared  to  A  =  9.5  pm, 
at  both  wavelengths  the  linear  relation  (7)  is  within 
a  factor  2.5  of  the  numerical  results  for  all  341  fog 
and  haze  size  distributions.  Examination  of  results 
at  numerous  other  infrared  wavelengths  show  that 
even  though  the  maximum  radius  condition  is  some¬ 
times  violated,  for  A  -  3.5-5. 3  pm,  8-10  pm  all 
numerical  results  for  the  341  fog  distributions  are 
within  a  factor  2.5  of  the  predictions  between 
absorption  and  liquid  water  content  given  by  rela¬ 
tion  (7)  and  listed  in  Table  1,  and  in  most  cases 
the  agreement  is  within  a  factor  2.  However,  at 
wavelengths  A  -  10.5- 12  pm,  the  numerical  results 
differ  by  as  much  as  a  factor  4  from  the  linear  (ela¬ 
tion  (7). 

Previously,  Platt  (197F)  found  an  approximate 
linear  relation  between  the  absorption  coefficient 

Tabu.  1.  At  a  given  wavelength  \  the  efficiency  facto  I  n 
absorption  (/„  (and  extinction  O,  I  can  K-  approximated  b-.  a 
stiaiehi  Ins  (J.,  «  '.<(<*,  ■  ex)  for  M/e  paiumetcrs  x  t  a  .  I  he 

Values  of  i.  and  o'  (ami  e)  are  determined  from  the  el,;  -jn.-y 
curves  (s-o.  eg.  1  igv  6  and  7|.  If  we  know  the  maximum  i.a-T  > 
r„  of  d'oplets  in  a  given  sbe  distribution,  the  table  gives  the 
wavelength  A  at  which  a  iinear  relationship  between  ahsoiption 
(or  extinsboni  and  liquid  water  con'ent  exists,  and  the  u— p: s>;. - i - 
at  va'ue  of  the  paran-e’er  / '  lor  c  l.  The  value  of  the  qu.ai: 
3-o'/':x/i  (or  3n<  -3v/il  multiplied  by  the  hon'd  water  coal.-nt  t- 
gives  t!:e  absorption  coolncient  o„  (or  exlinetien  coefficient  i 
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a,  at  A  -  II  /xm  and  liquid  water  content  IV'  of  non¬ 
precipitating  stratocumulus  clouds.  Platt  performed 
Mie  calculations  on  25  measured  cloud  droplet 
distributions,  plotted  the  values  of<r„  at  A  =  11  /xm 
vs  W,  did  a  least-squares  fit  through  the  resulting 
data  points  and  obtained  =  76.5  VV,  comparing 
modestly  well  with  our  Eq.  (7)  prediction  from  Table 
I  of  «r,  =  93  VV,  where  the  absorption  is  in  km-1 
and  the  liquid  water  content  in  g  m~3.  The  reason 
for  our  overprediction  of  the  absorption  is  that  a 
significant  number  of  droplets  in  Platt's  distributions 
have  radii  >7.2  /xm,  the  maximum  value  allowable 
in  the  Q„  =  c'x  approximation  at  A  =  11  /tm.  Had 
Platt  chosen  a  slightly  shorter  wavelength,  say, 
A  =  9.5  fxin,  he  would  have  found  an  even  better 
correlation  of  absorption  and  liquid  water  content, 
as  well  as  better  agreement  with  our  linear  predic¬ 
tion  (7),  since  the  maximum  radius  restriction  is 
more  nearly  satisfied  for  clouds  at  the  shorter 
wavelength. 

5.  Application  of  extinction — absorption-liquid  water 
content  relationships 

The  unique,  linear,  size-distribution-independent 
relationship  between  extinction  at  A  ~  11  /xm  and 
liquid  water  content,  and  between  absorption  at 
A  =  3.8,  9.5  gm  and  liquid  water  content  in  fog 
has  several  practical  '.pp’.ications.  For  example, 
for  a  fog  with  liquid  water  content  IV  =  0.1  g  m~3 


Fig.  8.  Variation  of  al  .wiPtion  coefficient  with  liquid 
water  content  in  atmospheric  fo  .  and  h.i/c  lor  341  size  distnbu- 
lion  measuiements  made  at  tiil’cicnt  googr.iphic  locales  and 
under  a  variety  of  mhM  cor  plop*  t.:*l  conditions.  In  the  infrared 
region  around  X  -  3  8  /urn  then  ousts  a  line.:.,  si/e-distribu* 
Cion- in  dependent  relator  the  vohmt.  ;iKorp:ion 

coefficient  ir,.  and  the  liquid  w.tfct  content  VI  oi  th*  form  of 
t.q  (7).  Consequently,  the  riMt'is  of  all  measuiements  arc 
close  to  a  straight  line.  I  he  pu-*!icted  relation  between  ab¬ 
sorption  (f,,  and  liquid  Wait:!  con'eni  W  accoiding  to  l.q  (7)  is 
shown  b>  the  straight  line 


Fig.  9.  As  in  Fig.  8  except  for  X  -  9.5  /em. 


we  predict  extinction  and  absorption  coefficients  of 
o>( A  =  11  gm)  =  13  km"1,  cr„(3.8  /tm)  =  1.8  km"1, 
cr„(5.3  /am)  =  3.1  km"1  and  o-u(9.5  /am)  =  5.7  km"1. 
Integrated  liquid  water  content  along  a  path  in 
fog  and  haze  could  be  inferred  from  a  CO.  laser 
(A  =  10.6  /xm)  transmissometer  measurement  ac¬ 
cording  to  (1).  Of  course  the  path  must  be  short 
enough  that  multiple-scattering  effects  and  forward¬ 
scattering  -orrections  (Deepak  and  Box,  1978V 
arc  negligible.  Carlon  ci  at.  (1977)  previously  realized 
this  applica.ion  for  a  transmission  measurement  at 
A  =  12.5  /xm.  An  application  of  the  crn-VV  rela¬ 
tionship  is  that  liquid  water  content  of  fog  at  a 
particular  point  could  be  determined  by  measure¬ 
ment  of  fog  droplet  absorption  in  the  spectral 
region  A  ~  3.8,  9.5  /xm  according  to  Eq.  (7).  The 
spectrophone  technique  has  been  demonstrated  to 
be  suitable  for  in  < si;u  measurement  of  particulate 
absorption  by  Bruce  and  Pinnick  (1977),  so  that 
measurement  of  fog  droplet  absorption  with  either  a 
CO,  laser  spectrophone  tuned  to  a  wavelength 
A  *=  9.5  /xm,  or  a  OF  laser  spectrophone  (A  ~  3.S 
/xm)  couiu  be  used  to  infer  fog  liquid  water  content. 
Of  course  the  spectrophone  measurement  of  ab¬ 
sorption  could  also  be  used  to  infer  fog  extinction 
at  A  ~  II  /urn.  1  he  Eq.  (7)  relation  between  fog 
absorption  and  liquid  water  content  might  also  be 
used  in  calculation  of  fog  emissivitics. 

6.  Conclusions 

Chylek's  (ll>7,S;  prediction  of  a  linear  relation, 
independent  of  the  form  of  the  size  distribution, 
between  exonetion  at  A  —  II  //in  and  liquid  water 
content  of  log  has  been  veiified  within  a  faetoi 
for  341  ditlcrcnt  tog  and  haze  droplet  disliihutiot 
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measured  under  a  variety  of  meteorological  condi¬ 
tions.  The  prediction  generally  works  better  for 
radiation  fogs  than  advection  fogs.  A  similar  linear 
relation  between  fog  droplet  absorption  at  A  ~  3.8, 
9.5  fim  and  liquid  water  content  has  been  derived 
and  validated  using  the  same  341  distributions. 
However,  there  exists  no  size-distribution-inde¬ 
pendent  relation  between  extinction  in  the  visible 
(A  =  0.55  few)  and  fog  liquid  water  content.  Three 
practical  applications  o?  these  findings  are  sug¬ 
gested:  1)  inference  of  fog-integrated  liquid  water 
content  along  a  path  by  measurement  of  laser 
transmission  (at  X  =*•  11  pm)  along  that  path;  2) 
inference  of  fog  liquid  water  content  at  a  particular 
point  from  measurement  of  fog  droplet  absorption 
with  a  DF  or  CO*  laser  spectrophone;  and  3) 
calculation  of  fog  emissivities  in  the  infrared  from 
knowledge  of  fog  liquid  water  content. 
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ABSTRACT 

Applying  perturbation  theory  we  have  derived  the  first  order 
perturbation  corrections  to  the  scattering  characteristics  for  the 
case  of  light  scattering  by  slightly  deformed  spheres.  Numerical 
results  are  compared  with  those  obtained  using  the  extended  boundary 
condition  method.  The  range  of  applicability  of  the  first  order 
perturbation  corrections  are  discussed. 


The  application  of  perturbation  theory  to  light  scattering  by 
non-spherical  particles  has  experienced  very  little  development  since 
it  was  first  derived  by  Yeh  (1964,  1965)  and  Erma  (1968).  It  is  the 
purpose  of  this  report  co  investigate  the  application  of  first  order 
perturbation  theory  to  light  scattering  by  various  non-spherical 
particles.  The  results  of  perturbation  theory  are  then  compared  with 
calculations  for  the  same  particles  using  the  extended  boundary  con¬ 
dition  method  (KLCM)  developed  by  Waterman  (1971),  and  Barber  and 
Yeh  (1075;. 


A  general  formalism  for  calculating  the  nt*1  order  perturbation 
term  has  been  developed  by  Yeh  (1964,  1965)  and  Erma  (1968).  The 
particle  is  assumed  to  be  of  the  form: 

r  =  rg  (l  +  e  f(6,o))  (1) 

where  rs  is  the  radius  of  an  unperturbed  sphere,  f(0,<f>)  describes  the 
shape  of  the  irregularity  on  the  sphere.  It  is  further  assumed  that 
e  «  1  and  that  |cf(0,^)]  <  1.  The  calculations  reported  herein 
used  the  following  functional  forms  for  f(0,<p): 


f(0,4>)  =  T2  (cos&) 

(2) 

f(e,4>)  =  T4(cos0) 

(3) 

where  Tn(cos0)  is  the  n*-^1  order  Chebyshev  polynomial. 

The  shapes  of  these  particles  are  illustrated  in  Figure  1.  The 
particle  is  assumed  to  be  oriented  such  that  it  is  rotationally 
symmetric  about  the  z-axis.  The  incident  light  is  assumed  to  be 
traveling  in  the  indirection  with  E  vector  in  the  x-direction.  First 
order  perturbation  theory  assumes  that  the  scattering  coefficients 
for  the  non- spherical  particle  can  be  expressed  as 


(4) 

(5) 


where  an,  bn  are  the  Mie  scattering  coefficients  for  a  sphere  of 
radius  rq;  a.-,,  b^  are  the  first  order  corrections  due  to  the  non- 
sphericity  of  the  particle. 

Calculations  were  performed  for  the  case  where  the  index  of  re¬ 
fraction  n  -  1.5  and  e  =  ±0.1.  Results  for  other  indexes  of  re¬ 
fraction  and  values  of  e  have  been  reported  in  Chylok  c t  at.  (1978). 
The  efficiencies  arc  calculated  by  normalizing  the  extinction  cross 
sections  by 

Figure  2  shows  the  extinction  efficiency  for  both  perturbation 
theory  and  K!'.C“  tor  the  case  of  f(0,^)  ~  T?(0)  and  c.  =  +0.1.  The 
extinction  efficiency  for  a  sphere  of  radius  rg  is  also  shown  in 
Figure  2.  ■ r  is  evident  from  Figure  2  that  perturbation  theory 

agrees  quite  will  with  1'BCM  for  Xj.  <  6.  For  values  of  xg  >  6  there 


FIG.  1.  The  two  shapes  considered  are:  a)  r  =  rg(l  +  e  T2(cos0)), 
b)  r  =  rs(l  ±  e  Ti,(cos0)). 


FIG.  2.  Tlie  extinction  effi¬ 
ciency  for  light  scattering  by 
a  particle  whose  shape  is  given 
by  r  =  rs(l  +  0.1  T?(cosO)J. 

The  solid  line.  ( - )  is  for  a 

sphere  with,  r  =  rs.  The  dashed 

line  ( - )  represents  the  exact 

results  of  the  EBCM.  The  dotted 
line  (•••)  represents  the  re¬ 
sults  of  first  order  perturba¬ 
tion  theory. 


is  very  poor  agreement  between 
perturbation  theory  and  EBCM  due 
to  the  narrowness  of  the  resonance 
in  the  extinction  efficiency.  The 
effect  of  the  width  of  a  resonance 
on  perturbation  theory  has  been 
discussed  by  Chylek  et  al.  (1978). 
There  it  was  shown  that  in  the 
resonance  region  one  would  expect 
perturbation  theory  to  apply  when 

|  e |  <  r/2  xg  ,  (6) 


where  T  is  the  half  width  of  a 
Lorentzian  shaped  resonance. 
Equation  (6)  implies  that  for 
large  values  of  Xj.,  where  T  is 
expected  to  be  small,  one  would 
need  a  very  small  e  for  perturba¬ 
tion  theory  to  agree  with  the  exact 
calculations. 

Figure  3  shows  the  extinction 
efficiency  for  the  same  shape 
function  but  with  e  =  -0.1.  The 
observations  made  above  for  Figure 
2  apply  to  this  case. 


f(0,*)  =  T2(0>  '\wr 

rn  =  1.5  €  =  +0.1 
. Perturbotion  Theory 


The  difference  between  perturbation  theory  and  EBCM  can  best  be 
illustrated  by  considering  the  relative  difference  between  the  re¬ 
sults  obtained  using  these  two  methods.  This  difference  is  defined 

as 


Pert  _  EBCM 
^ext  “  ^ext 


Q 


EBCM 

ext 


(7) 


which  is  a  function  of  xg.  Figure  4  shows  the  absolute  value  of 
equation  (7)  for  the  case  of  e  =  +0.1  and  f  (0 ,4>)  =  12(6).  Since 
terms  of  order  e2  have  been  neglected  in  these  calculations  an  error 
of  10~2  would  be  expected.  It  can  be  seen  from  Figure  4  that  the 
difference  lies  between  10-3  and  10-1.  A  striking  feature  of  this 


FIG.  3.  The  extinction  effi¬ 
ciency  for  light  scattering  by 
a  particle  whose  shape  is  given 
by  r  =  rs(l  -  0.1  T2(cos9)J. 

The  solid  line  (— — )  is  for  a 
sphere  of  radius  r  =  rs.  The 

dashed  line  ( - )  represents 

the  exact  results  of  the  EBCM. 

The  dotted  line  (•••)  repre-  FIG.  4.  The  absolute  value  of 

sents  the  results  of  first  the  relative  difference  between 

order  perturbation  theory.  perturbation  theory  and  the  EBCM 

for  the  particle:  r  =  rs(l  + 

0.1  T2CCOSO))  as  a  function  of 
xg  =  2irrs/X. 


figure  is  the  oscillatory  behavior  of  |ej.  These  oscillations  are 
due  to  the  resonance  effect  discussed  above  for  Figures  2  and  3. 

There  exists  a  one-to-one  correspondence  between  the  position  of  the 
maxima  in  relative  difference  in  Figure  4  and  the  position  of  the 
resonance,  of  the  extinction  curve.  An  overall  growth  in  absolute 
relative  difference  can  be.  observed  for  increasing  values  of  xs  due 
to  the  narrowing  of  the  resonances. 

Similar  calculations  for  f(0,<j>)  =  1'4(0)  and  c  =  +0.1  are  shown 
in  Figure  5.  For  this  case  little  difference  can  be  seen  to  exist 
for  xg  <  7.  It  is  also  observed  that  the  T[,(0)  particle  has  a 
smaller  extinction  than  the  sphere  for  this  size  range.  Figure  6 
illustrates  the  extinction  efficiency  for  a  shape  function  f(0,<)))  = 
Tj^fO)  and  e  =  -0.1. 

Figure  7  shows  the  absolute  value  of  the  relative  difference  for 
the  particle  with  e  ■=  +0.1  as  a  function  of  x^.  The  relative 

difference  from  this  figure  lies  between  2.5  x  10-1*  to  2.2  x  10"'. 

The  oscillatory  behavior  of  this  difference,  correlates  with  the  res¬ 
onances  in  Figure  5. 


FIG.  5.  The  extinction  effi¬ 
ciency  for  light:  scattering  by 
a  particle  whose  shape  is  given 
by  r  =  rs(l  +  0.1  Tt  (eosfl)J  . 

The  solid  line  ( - )  in  for  a 

sphere  of  radius  r  -  r,_.  The 

dashed  line  ( - )  represents 

the  exact  results  of  the  ERCM. 
The  (letted  line  (•••)  represents 
the  results  of  first  order  per¬ 
turbation  theory. 


FIG.  6.  The  extinction  efficiency 
for  light  scattering  by  a  particle 
whose  shape  is  given  by  r  -  rs 
(1  -  0.1  T  i(  (cos  0)  .  The  solid 
line  (• — )  is  for  a  sphere  of  radius 

r  -  r;;.  The  dashed  line  ( - ) 

represents  the  exact  results  of 
the  K.VM.  The  dotted  line  (•■•) 
represeu  s  the  results  of  first 
order  perturbation  booty 


In  conclusion,  this  study 
indicates  that  first  order  per¬ 
turbation  theory  as  applied  to 
light  scattering  by  non-spherical 
particles  is  limited  to  small  de¬ 
formations  and  low  values  of  the 
size  parameter.  The  limitation 
arises  because  of  the  breakdown 
of  first  order  perturbation  theory 
in  the  region  of  resonances  of 
the  scattering  coefficients. 
However,  these  problems  can  be 
eliminated  by  taking  higher  order 
terms  in  the  perturbation  ex¬ 
pansion.  The  computational  feasi¬ 
bility  of  such  calculations  is 
possible  since  these  calculations 
require  little  storage  and  com¬ 
puter  time. 


FIG.  7.  The  absolute  value  of 
the  relative  difference  between 
perturbation  theory  and  the  EBCM 
for  the  particle:  r  =  rs(l  +0.1 
T4(cos0)J  as  a  function  of  xg  = 
2itrs/X. 
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ABSTRACT 

Resonances  in  partial  wave  scattering  amplitudes  an  and  bn  are 
responsible  for  the  ripple  structure  ot  the  extinction  curve  and  for 
sharp  peaks  in  the  backscattering  (glory) .  A  connection  between  the 
resonances  and  surface  waves  is  suggested. 


1.  RESONANCES  IN  PARTIAL  WAVE  AMPLITUDES  an  AND  bn 

It  is  convenient  to  write  the  Hie  partial  wave  amplitudes  an 
and  bn  in  the  form 


a 

n 


A 

n 


<ls>  *  b„  '  <lb> 

n  n 


where  7n,  C,-,,  Bn  and  D.,  are  real  functions  for  real  refractive  index 
m.  Their  explicit,  form  has  been  given  elsewhere  (Chylek,  1973). 


For  simplicity  of  discussion  we  consider  a  real  index  of  refraction 
m.  Then  it  follows  that  the  Re{an)  and  Re{bn}  reach  a  series  of 
absolute  maxima  Re{an)  =  1  and/or  Re{bn}  =  1  at  such  values  of 
x  (x  =  2irr/X  is  a  size  parameter)  where  In{anJ  =  0  and/or  Im{bn}  =  0. 
At  sufficiently  high  values  of  n,  the  first  maximum  in  each  partial 
wave  amplitude  nn  and  bn  has  the  form  of  a  sharp  peak.  In  Figures 
1A  and  IB  the  first  peaks  in  Re{ajQ)  and  Refa^}  are  shown,  and 
Figures  2A  and  2B  show  the  first  sharp  peak  in  Re{b3)}  as  well  as 
Im{b3i).  These  sharp  peaks  we  call  resonances.  Wo  notice  that  with 
increasing  n  the  first  peak  in  each  partial  wave  becomes  narrower; 
also  with  an  increasing  value  of  refractive  index  m  the  peaks  become 
narrower.  At  higher  values  of  n  not  only  the  first  peak  but  also  a 
few  of  the  following  peaks  may  be  sharp  enough  to  be  called  reson¬ 
ances. 


2.  WHAT  ARE  SURFACE  WAVES? 

There  seems  to  be  no  generally  accepted  definition  of  what  is  a 
surface  wave  in  the  case  of  scattering  of  electromagnetic  waves  on 
small  particles.  We  are  going  to  define  what  we  mean  by  a  surface 
wave.  We  have  no  intention  in  claiming  that  our  definition  should  be 
generally  accepted.  We  simply  feel  that  since  we  deal  with  surface 
waves,  we  have  to  define  clearly  what  is  meant  by  it  in  this  paper. 

We  make  a  reference  to  the  two  well-known  phenomena  generally 
attributed  to  surface  waves,  namely  the  ripple  structure  of  the  ex¬ 
tinction  cross  section  (Fig.  3A)  and  to  the  glory  phenomena  of  back- 
scattering  (Fig.  3B).  From  a  large  number  of  published  papers  and 
from  available  monographs  on  light  scattering  we  can  learn  that  the 
ripple  structure  in  the  extinction  cross  section  is  presumably  caused 
by  the  interference  phenomena  between  a  forward  diffracted  wave  and 
the  surface  wave.  Similarly  the  peaks  appearing  in  the  backscattering 
cross  section  (glory)  are  supposed  to  be  a  result  of  interference 
between  the  backward  reflected  wave  and  the  surface  wave.  Without 
trying  to  discuss  the  validity  of  the  above  statements,  we  will  show 
that  both  the  ripple  structure  of  extinction  cross  section  and  the 
glory  in  the  backscattering  are  caused  by  the  previously  described 
resonances  in  the  partial  wave  amplitudes  an  and  bn. 

Let  us  go  back  to  Figures  2A,  3A  and  3B.  We  notice  that  the 
resonance  ii:  b 3 j  (Fig.  2A)  occurs  somewhere  between  26.9  <  x  <  27.0. 

We  also  notice  that  in  the  same  region  of  x  there  appears  a  sharp 
peak  in  the  ripple  structure  (Fig.  3A)  and  a  glory  (sharp  peak)  in 
the  backscattering  cross  section  (Fig.  3B) .  This  leads  to  a  con¬ 
jecture,  suggested  some  time  ago  by  one  of  the  authors,  that  there 
is  one-to-one  correspondence  between  the  partial  wave  resonances  and 
the  ripple  structure  in  the  extinction  and  the  glory  in  the  back- 
scattering. 


3/ 


FIG.  1A.  The  first  resonance  in  the  partial  wave  ajg  occurs  at 
x  =  9.203  f,or  refractive  index  m  =  1.50. 


FIG.  IB.  Vith  increasing  n  the  resonances  become  narrower.  The 
first  resonance  in  a?g  occurs  at  x  =  16.65. 


FIG.  2A.  In  general  the  an  resonances  are  sharper  than  the  bn  res¬ 
onances  of  the  same  value  of  n.  Also,  with  higher  refractive  index 
m,  resonances  are  narrower.  The  first  resonance  in  b3i  for  refrac¬ 
tive  index  m  =  1.33  has  a  width  comparable  to  the  width  of  a2o 
(with  m  -  1 .50) . 


FIG.  2B.  The  imaginary  part  of  b^j  goes  through  a  zero  in  the  di¬ 
rection  f rou  positive  to  negative  values  at  the  x  value  at  which  the 
real  part  of  b 3 2  readies  its  maximum  value  Re{b3j}  =  1. 


FIG.  3A.  The  normalized  extinction  cross  section  QgxT  (efficiency 
for  extinction)  shows  a  ripple  structure.  Each  peak  in  a  ripple 
corresponds  to  a  definite  resonance  in  the  partial  waves  an  or  bn 
(surface  waves).  The  peak  between  26.9  <  x  <  27.0  corresponds  to 
the  631  resonance. 


FIG.  3B.  The  backscaitcring  is  completely  dominated  by  the  surface 
waves  (partial  wave  resonances  in  an  and  bn) .  Again  the  peak  in  the 
backscattering  between  26.9  <  x  <  27.0  is  caused  by  he  first  reson¬ 
ance  in  the  b 3 1  partial  wave. 
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FIG.  4A.  To  demonstrate  that  the  ripple*  structure'  of  the  extinction 
curve  is  caused  by  the*  partial  wave  resonances  in  a  and  bn,  we  have 
calculated  a-rvia  (ijyp ,  however,  now  we  have  set  l.> 3 1  -  0  in  the  reson¬ 
ance  region  24.6  <  x  <  27.4.  Comparing  Figs.  3A  and  4A  we  notice  that 
the  peak  in  Q.-v-p  between  26.9  <  x  <  27.0  disappears,  which  shows 
clearly  that  the-  peak  is  formed  by  the  l> 3 j  resonance. 


FIG.  4H.  When  the  b.j  resonance  is  n  „>oved  f  1  on  tin*  backseat  t  or  i  up, 
colonial  icr  the  peak  in  t  he  backseat l or in>  (glory)  between  36.9  <  x  «■ 
27.0  disap;'-,  rs,  which  shows  that  the  peak  is.  formed  by  the  b,j  1  osi  1 
fine  e . 


^jgj 


Let  us  consider  the  resonance  contributions  of  3)q(x)  and  a2o(x)- 
The  peak  of  the  extinction  curve  at  x  =  9.203  (Kip,.  5A)  corresponds 
to  the  ajQ  resonance  and  the  peak  at  x  =  26.65  (Fig.  6A)  corresponds 
to  the  a?o  resonance.  Using  Equation  (4)  we  can  calculate  the  back¬ 
ground  contributions  (see  Fig.  5B  and  6B)  and  we  obtain  Qt.-\gj>  }> 

(9.203)  =  1.9U.'1,  Qj.-v-j-  j» ( 1 6 . 65)  -  2.502.  Subtracting  the  values  from 
the  extinctions,  Qf;xT *  (9 . 203 )  =  2.475  and  QgX'f  (16.65)  =  2.798,  we 
obtain  for  the.  resonance  contributions  Q]  \;t  RES (9 . 203)  =  0.493  and 
^EXT  RKS^^  *  66)  ~  296  or,  in  the  form  of  a  ratio  of  the  a2Q  and 

ajo  resonance  contributions,  we  obtain 


^ EXT , RES  ^16,65^ 

<Wes  <9-203) 


0.600  . 


If  we  calcul.  t ■  the  ratio  of  the  same  resonance  contributions  using 
an  approximate  form  of  equation  (3)  we  obtain 


,_KES 
,  RES 


(16.65) 

(97203) 


41  (9.203)2 
21  (16.65F 


0.596 


in  excel  1  en*  . . .  with  the  previously  obtained  value  of  0.600. 


We  coco  that  the  resonance  or  surface  wave  contribution  to 
the  ext  in'. !  described  by  equation  (3)  is  in  agreement  with 

numeric:.!  r. 


4.  SURFACE  WAVES  IN  EXPERIMENTAL  MEASUREMENTS 

Recent  1 >  a  technique  of  optical  levitation  was  used  to  support 
a  liquid  droplet  in  a  laser  beam  (Ashkin,  1970;  Ashkin  and  Dzicdzic, 
3976,  19//).  One  of  the  authors  (A.  Ashkin)  used  this  technique  in 
connect  io. i  .  .  a  a  tunabl e  laser  to  measure  the  hackscattering.  An 
example  ot  i.  ■  a  sure  ment  of  the  backseat  l  e  red  intensity  is  shown 

in  Figure  Three  strong  peaks  observed  in  the  wavelength  interval 

correspond  :  >  :  io  glory  effect.  Since  the  droplet  radius  is  known 
only  appro.  ■  it  is  not  possible  to  label  the  horizontal  axis 

in  size  p;  r  visits  at  this  time. 

The  f  .i  .  .  p  'nt  v.v  want  to  demonstrate  is  that  there  is  only 
one  way  to  id  ;  f ...  the  experimental  measurement  with  theoretical 
cn1cel.it  :  •..hill  tiie  droplet  radius  is  essentially  unknown. 

The  rests  ,  :  t  be  backseat  ter  i  ng  calcul.:'.  ion::  with  t  he  size  para¬ 
meter  inc:  •  .St  if  i\.\  -  0.001  are  shown  in  Figui  e  7b.  We  notice  how 
the  sh.iri'  r  •  •>.•  e  (glory)  is  super  imp.-,  e.l  on  a  broader  background 
(really  <  ,  ae  ba.  kg, round  con  ing,  fror:  the  teims  with  u  <  x  plus 


FIG.  5A.  The  normalized  extinction  cross  section  Qg^x  is  composed 
from  the  resonance  contributions  superimposed  on  the  background  term. 


FIG.  5B.  When  the  first  ajo  resonance  is  removed  around  x  =  9.2 
only  the  background  term  remains. 


FIG.  6A.  The  peak  in  Qgxj  around  x  =  16.6  is  caused  by  the  a2o 
resonance. 


FIG.  6B.  When  the  ajo  resonance  is  removed  from  the  calculations 
(we  set  a?£  =  0)  the  peak  around  x  =  16.6  disappears  and  only  the 
background  term  is  left. 


Ut: 
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FIGS.  7A-/G.  The  experimental  measurement  (A)  of  the  backseat tering 
shows  resonances  superimposed  upon  the  background.  We  notice  narrow 
resonance  peaks  as  well  as  the  higher  order-  broader  resonances.  The 
distance  between  the  neighboring  resonances  of  the  same  order  depends 
slightly  on  the  order  of  resonances  and  on  the  size  parameter.  Be¬ 
cause  of  chat ,  each  group  of  resonances  (denoted  by  1,  2  and  3)  has 
a  different  characteristic  structure,  i.c.,  different  morphology. 

The  same  structure  appears  in  a  coarse  numerical  calculation  (B) , 
and  it  is  just  this  characteristic  morphological  structure  of  each 
group  of  resonances  which  allows  us  unambiguously  to  identify  the  ex¬ 
perimentally  observed  peaks.  A  detailed  numerical  calculation  (C) 
reveals  additional  narrow  resonances  that  have  not  been  observed  ex¬ 
perimentally  because  their  width  is  below  the  resolution  power  of  the 
experimental  equipment,  used.  The  refractive  index  of  oil  droplet 
used  is  m  1.42. 


a  broader  higher  order  resonances  in  partial  waves  with  n  =  x).  We 
conclude  that  the  morphology  of  the  resonance  structure  leads  to  an 
unambigous  identification  of  the  experimentally  measured  peak  1  in 
Figure  7  A  with  the  theoretically  calculated  peak  1  in  Figure  7:  . 
Now,  since  we  know  precisely  the  wavelength  >.  at  which  peak  1  was 
observed  as  well  as  the  size  parameter  x  at  which  peak  1  occurs  in 
the  numerical  calculations,  we  can  determine  with  a  high  degree  of 
accuracy  the  radius  of  the  observed  droplet,  as  has  been  already 
demonstrated  elsewhere  (Chylek  et  al.  ,  1978a,  1978b). 


We  were  fully  aware  that  by  using  the  step  size  of  Ax  =  0.001 
in  our  numerical  calculations  we  could  lose  many  of  the  narrow  peaks 
in  the  backscattering  structure.  By  using  a  computer  resolution  of 
the  same  order  as  the  resolution  of  the  experimental  equipment  used, 
we  hoped  to  get  comparable,  structures  as  shown  in  Figures  7A  and  7B. 
By  increasing  the  resolution  of  our  numerical  calculations  we  obtain 
the  full  backscattering  structure  (Fig.  7 C).  It  follows  that  the 
resolution  power  of  the  instrumental  arrangement  used  was  not  high 
enough  to  observe  the  first  and  the  second  order  resonances.  To 
observe  those  narrow  resonances  the  resolution  power  cf  experimental, 
observation  has  to  be  increased  by  several  orders  of  magnitude. 


5.  SURFACE  WAVES  ON  NON-SPHEKICAL  PARTICLES 


Our  preliminary  results  indicate  that  surface  waves  do  exist  at 
least  in  definite  orientations  of  axially  symmetric,  ron-spherical 
particles.  On  the  other  hand,  the  strength  of  surface  waves  is  con¬ 
siderably  reduced  when  non-sphericai  particles  are  in  random  orien¬ 
tation.  The  jo  conclusions  were  deduced  from  numerical  calculations 
using  the  extended  boundary  condition  method.  Detail;,  of  these  cal¬ 
culations  wiki  be  reported  elsewhere. 
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